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Abstract. Continuous variable remote state preparation and teleportation are 
analyzed using Wigner functions in phase space. We suggest a remote squeezed 
state preparation scheme between two parties sharing an entangled twin beam, 
where homodyne detection on one beam is used as a conditional source of 
squeezing for the other beam. The scheme works also with noisy measurements, 
and provide squeezing if the homodyne quantum efficiency is larger than 50%. 
Phase space approach is shown to provide a convenient framework to describe 
teleportation as a generalized conditional measurement, and to evaluate relevant 
degrading effects, such the finite amount of entanglement, the losses along the 
line, and the nonunit quantum efficiency at the sender location. 



1. Introduction 

Let us consider an entangled state described by a density matrix i? on a bipartite 
Hilbert space Hi (^0.2- A measurement performed on one subsystem reduces the 
otlier one according to the projection postulate. Each possible outcome, say x, occurs 
with probability , and corresponds to a different conditional state 

Px = Tri2 [R n, ® h] , & = — Tn [R ® h] . (f ) 

Px 

Tlx is the probability measure (POVM) of the measurement (acting on the Hilbert 
space of the first subsystem) and I2 the identity operator on the second Hilbert space. 
Tri2 [...] denotes full trace, whereas Tr^ [...], j = 1,2 denotes partial traces. 

Eq. Q shows that entanglement and conditional measurements can be powerful 
resources to realize (probabilistically) nonlinear dynamics that otherwise would not 
have been achievable through Hamiltonian evolution in realistic media. Since 
entanglement may be shared between two distant users (the sender performing 
the measurement, and the receiver observing the conditional output), the inherent 
nonlocality of entangled states permits the remote preparation of the conditional states 
Qx, a protocol that may be used to exchange quantum information between the two 
parties sending only classical bits A different kind of remote state preparation 
is teleportation [2], where the measurement depends on an unknown reference state 
which may be recovered at the receiver location independently on the outcome of the 
measurement. 
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In this paper, we focus our attention on continuous variable (CV) remote state 
preparation . In particular, we analyze in detail an optical scheme for remote 
preparation of squeezed states by realistic (noisy) conditional homodyining. Our 
analysis is based on a phase-space approach, and this is motivated by the following 
reasons: i) entanglement in optical CV quantum information processing is provided by 
the so-called twin-beam (TWB) state of two field modes |A)) = \/l — J2p \p)\p)^ 
< A < 1; the corresponding Wigner function is Gaussian; ii) trace operation 
corresponds to overlap integral [3, and the Wigner function of (realistic) homodyne 
POVM is also a Gaussian. By Wigner calculus we will be able to derive simple 
analytical formulas for conditional outputs, also in the case of noisy measurement 
at the sender location. In addition, we will show that phase-space approach is a 
convenient framework to describe CV teleportation as a conditional measurement , 
and to evaluate relevant degrading effects, such the finite amount of entanglement, 
the losses along the transmission channel, and the nonunit quantum efhciency at the 
sender location. 



2. Conditional measurement in phase space 



TWB is the maximally entangled state (for a given, finite, value of energy) of 
two modes of radiation. It can be produced either by mixing two single-mode 
squeezed vacuum (with orthogonal squeezing phases) in a balanced beam splitter [3] 
or, from the vacuum, by spontaneous downconversion in a nondegenerate parametric 
optical amplifier (NOPA) The evolution operator of the NOPA reads as follows 
Ur = exp [r (a^6^ — a6)] where the "gain" r is proportional to the interaction-time, 
the nonlinear susceptibility, and the pump intensity. We have A = tanhr, whereas 
the number of photons of TWB is given by N = 2sinhV = 2AV(1- A^). In view of 
the duality squeezing/entanglement via balanced beam-splitter 6 the parameter r is 
sometimes referred to as the squeezing parameter of the twin-beam. Throughout the 
paper we will refer to mode a as "mode 1" and to mode b as "mode 2". The Wigner 
function VF[twb](xi, a;2, ?/2) of a TWB is Gaussian, and is given by (we omit the 
argument) 



exp 



(Xi+X2f {Vl+y2f {Xl-X2f (2/1-2/2) 



21 



4a2 



where the variances are given by 

2 1 
0-+ = 4exp{2r} 



4a? 



1 

-exp{-2r} 



4a? 



(2) 



(3) 



Specializing Eq. for i? = |A))((A| we have 

p, - ((A|n,§5/2|A» - (l-A2)Tri 

= — Tri[|A))((A|n,®/2] , 

Px 

where, in the expression of Px , we have already performed the trace over the Hilbert 
space 7^2- In the following, the partial traces in Eq. Q will be evaluated as overlap 
integrals in the phase space. The Wigner function of a generic operator O is defined 
as the following complex Fourier transform 

J2, 



W[0]{a) = j i2e"7-a7 Tr [O D{-i)] , 



(4) 
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where a is a complex number, and D('-f) — e'^°'^~^°' is the displacement operator. The 
inverse transformation reads as follows [7] 

= j (fa W[0]{a) e-2|«l' e^""^' {-f" e^"'' (5) 

Using the Wigner function the trace between two operators can be written as 

Tr[Oi O2] = TT / d'pWpim W[02m . (6) 



2.1. Remote squeezed states preparation 

Let us consider the optical scheme depicted in Fig. ^ A TWB is produced by 
spontaneous downconversion in a NOPA, and then homodyne detection is performed 
on one of the two modes, say mode 1. The POVM of the measurement, assuming 
perfect detection i.e. unit quantum efficiency, is given by 

\x)'s being eigenstates of the quadrature operator x — l/2(a + a^). The Wigner 
function of the POVM is a delta function 

W[U^]{xi) ^ S{xi - x) , (8) 

whereas that of the term A''^" in the first of Eqs. is given by 

{l-X')W[X^''^]{x,,y,) = {2naY'expl^-^^±^^ , (9) 

where the variance a depends on the number of photons of the TWB — + N). 
Using Eqs. (|HJl and JHJl it is straightforward to evaluate the probability distribution 

Px = JJJJdxidyidx2dy2 W[TWB]{xi,yi;x2,y2)W[U^]{xi) 
= (1 - A2)y"y" dxidyi W[X'''''](xi,yi) W[U^]{xi) 



= (2vra^)-^/^exp|-^| , (10) 
and the Wigner function of the conditional output state 

W[g^]{x2,y2) = dxidyi W[TWB]{xi,yi; X2,y2)W[Il^]{xi) 



(2.I]^2.S^)-V2exp^-I:£E_p)!_ . (u) 



The parameters in Eq. are given by 



= — x Li = Zj, = -(1 + A*) . (12) 

Eqs. pi|) and (|12|) say that Qx is a squeezed-coherent minimum uncertainty state of 
the form Qx — D{ax)S{rx)\G) i.e. a state squeezed in the direction of the measured 
quadrature Aa;^ = 1/46"^'"^, with squeezing parameter given by r^ = l/21og(l + N). 
Notice that this result is valid for any quadrature and therefore 
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the present scheme, by tuning the phase of the local oscillator in the homodyne 
detection, is suitable for the remote preparation of squeezed states with any desired 
phase of squeezing. Of course, we have squeezing for if and only if TV > i.e. if 
and only if entanglement is present. 

A question arises whether or not the remote preparation of squeezing is 
possible with realistic homodyne detection, i.e. with noisy measurement of the field 
quadrature. The POVM of a homodyne detector with quantum efficiency 77 is a 
Gaussian convolution of the ideal POVM 

n., = / exp{-fc|l!} \y){y\ , (13) 

with (7^ = \{i — r])/v ^ ■ The corresponding Wigner function is given by 

W[Il„,]{x,) = (27ra2)-V2exp|-i^i_^| . (14) 

Using (|14|l one evaluates the probability distribution and the Wigner function of the 
conditional output state, one has 

„2 



= [271(^2 +a2)]-i/2cxp^-^^-^^-5y|> (15) 



W[g.,]ix,,y,) ^ (2;rI]j„2vrI]^„)-V2exp|- (^^^^°-''^' - ^ } , (16) 
where 



WiV(jV + 2) _l l + jV(l-r;) ^2 

As a matter of fact, the conditional output g^ri is no longer a minimum uncertainty 
state. However, for 77 large enough, it still shows squeezing in the direction 
individuated by the measured quadrature i.e. Ax^ < 1/4. In order to obtain the 
explicit form of the conditional output state from the Wigner function W[gxi-i]{x2,y2) 
of Eq. ifT^ we use Eq. ^ arriving at 

= D{ax,,)S{r 

XT] )^th {r.,n)D\axn) , (18) 

where Vth = (1 + n-th)^^ EpI^-fi/ll + ^th)]^ \p){p\ is a thermal state with average 
number of photons given by 




l + rjN 

and the squeezing parameter is given by 

(20) 

1 (l + jV)(i + ^iV) 
'^^" = 4^°^ i + ■ (21) 

We have squeezing in gx-q if '^ijj < 1/4, and this happens for r] > 50% independently 
on the actual value x of the homodyne outcome. The values of efficiency that can 
be currently realized in a quantum optical lab is far above the 50% limit, and thus 
we conclude that conditional homodyning on TWB is a robust scheme for the remote 
preparation of squeezing. 
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2.2. Teleportation as a generalized conditional measurement 

The scheme for optical CV teleportatfon is depicted in Fig. |21 One part of a TWB is 
mixed with a given reference state cr in a balanced beam splitter, and two orthogonal 
quadratures x = are measured on the outgoing beams 

by means of two homodyne detectors with local oscillators phase-shifted by tt/2. The 
other part of the TWB is then displaced by an amount —a ~ —x — iy that depends on 
the outcome of the measurements, and the resulting state (averaged over the possible 
outcomes) is the teleported state. 

Overall, the measurement performed on the TWB is a generalized double 
homodyne detection [51 I10| (equivalent to generalized heterodyne), which can be 
described by the POVM [TOIIIII 

n„ L>(a)f7^L>^(a) , (22) 

...^ denoting transposition. Therefore, using Eq. (PJ, one has 



Pa 



((A|n„ ® /2|A)> = (1 - A2)Tri [a"^'" D{a)a^D\a) 



go. = —D{-a)TTi [|A))((A| D{a)a^ D^a) <g> h] D\-a) , (23) 

Pa 

while the teleported state is given by 

Jd^apa^Qa^ J d^ a D{~a)TTi [\\)) {{X\ D{a)a^ D^a) (g) h] D'' {-a) . (24) 

Using Wigner functions and taking into account that for any density matrix 

W[g^]{x,y) = W[g]ix,-y) (25) 
W[Dia)gD\a)]ix, y) = W[g]{x ~x^,y-y^), (26) 

with Xa = Re [a] and ya = Im[Q;], one has 

W[g]{x2,y2)^ 

= JJdxidyi JJdxadya W[TWB]{xi,yi- X2 + Xa,y2 + Ua) W[a-]{xi - Xa,-yi - Va) 

= JJdxidyi W[a]{xi,yi) JJdxadya W[twb]{xi + Xa, -yi - ya]X2 + Xa,y2 + Va) 
dxidyi r {xi -X2f + {yi -yif \ . 

exp(-^i±^l W[D{a,)aD\a,)]{x2,y2) , (27) 



TTKt. K 



with ai = xi + iyi and — exp{— 2r}. From Eqs. (|27|l and lO one has that the 
teleported state is given by 

g^ exp{-M!} D{a)<jD^a) , (28) 

which coincides with the input state only in the limit r — > oo i.e. for infinite energy 
of the TWB. Eq. H28|) shows that CV teleportation with finite amount of entanglement 
is equivalent to a thermalizing channel with thermal photons: this results has been 
obtained also with other methods However, the present Wigner approach may be 
more convenient in order to include other degrading effects such the nonunit quantum 
efficiency at the sender location and the losses along the transmission channel. 
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Nonunit quantum efficiency at the tiomodyne detectors affects ttie POVM of the 
sender, which become a Gaussian convoiution of the ideai POVM 

= / exp{- ' P } , (29) 

with — ~ rj)/rj 8 . On the other hand, fosses along the hue degrade the 
entanglement of the TWB supporting the tefeportation. The propagation of a TWB 
inside optical media can be modeled as the coupling of each part of the TWB with 
a non zero temperature reservoir. The dynamics can be described in terms of the 
two-mode Master equation 

^=Cgt= r(l + M)L\a]gt + T{1 + M)L\b]gt + rML\a'']gt + TML\b^gt (30) 
at 

where gt = git), T denotes the (equal) damping rate, M the number of background 
thermal photons, andLfO] is the Lindblad superoperator L[0]pt = OgtO^ — ^O^ Ogt — 
^gtOO^ . The terms proportional to L[a] and L[b] describe the losses, whereas the 
terms proportional to L[a^] and describe a linear phase-insensitive amplification 
process. This can be due either to optical media dynamics or to thermal hopping; 
in both cases no phase information is carried. Of course, the dissipative dynamics 
of the two channels are independent on each other. The master equation (|30|l can 
be transformed into a Fokker- Planck equation for the two-mode Wigner function of 
the TWB Using the differential representation of the superoperators in Eq. H3U|) the 
corresponding Fokker-Planck equation reads as follows 



drWr = 




Wr , (31) 



where r denotes the rescaled time t = (r/7) i, and 7 ~ 2M+1 drift term. The 
solution of Eq. (|31|l can be written as 

2 

Wr= Jdx[ JdxUdy[ Jdy', W[TWB]ix[,y[;x'„y',)f[Grix,\xr)GAyM) (^2) 

where VF[twb] is initial Wigner function of the TWB, and the Green functions 
Gr{xj\x'j) are given by 

(Xj - x' 6-57^)2 



GTixjlx'A = , exp 



D'^^il-e'-'n. (33) 



The Wigner function Wr can be obtained by the convolution (|32|l . which can be easily 
evaluated since the initial Wigner function is Gaussian. The form of Wr is the same 
of Ty[TWB] with the variances changed to 

al ^ (e-^V2 + D^) ^ (e-^V2 + D^) . (34) 

Inserting the Wigner functions of the blurred POVM IIq^ and of the evolved TWB 
in Eqs. (|23|l and H24() we obtain the teleported state in the general case, which is still 
given by Eq. (|28ll . with the parameter now given by 



2 > e-"-2r_^_(2M + !)(!_ e-rt) + A2 . (35) 

Eqs. (|28|l and (|35|l summarize all the degrading effects on the quality of the 
teleported state. In the special case of coherent state teleportation a = \z){z\ 
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(which corresponds to original optical CV teleportation experiments the fidelity 
F — {z\g\z) can be evaluated straightforwardly as the overlap of the Wigner functions. 
Since VF[z](q!) — 2/Tr e~^l"~^'^ is the Wigner function of a coherent state we have 



1 + e-2'-rt + (1 - e-r*)(2A/ + 1) + {1 - Tj)/rj ' 

The condition on the fidelity, in order to assure that the scheme is a truly nonlocal 
protocol, is given by F > 1/2 0, i.e. 

e-^r-Tt ^ _ g-rt)(2Af + 1) + (1 - 77)^ < 1 . 

Therefore, the bound on the quantum efficiency to demonstrate quantum teleportation 
is given by 

1 

^ 2 - 6-2-^-^* - (1 - e-r*)(2Af + 1) ' 
If the propagation induces low perturbation i.e. if F ~ and M ~ we have 
77 > (2 — e~2»')-i^ which ranges from 1/2 to 1, and represents the range of "useful" 
values for the quantum efficiency. If T and M are not negligible then, for the same 
initial squeezing, we need a larger value of the quantum efficiency. Moreover, since 
quantum efficiency should be lower or equal to unit 77 < 1, we may derive a bound on 
the initial squeezing that allows to demonstrate quantum teleportation. This reads 
as follows e~2'' < (2Af + 1) — 2Me^* . Remarkably, if the number of thermal photons 
is zero, i.e. if the TWB is propagating in a zero temperature environment, then any 
value of the initial squeezing parameter make teleportation possible, of course if the 
quantum efficiency at the receiver location satisfies ?? > (2 — e~2''-r< — 1 + e'"*)~^. 

3. Conclusions 

A method for the remote preparation of squeezed states by conditional homodyninig 
on a TWB has been suggested. The scheme has been studied using Wigner function, 
which is the most convenient approach to describe effects of nonunit quantum efficiency 
at homodyne detectors. The method is shown to provide remote squeezing if the 
quantum efficiency is larger than 50%. Since downconversion correlates pair of modes 
at any frequencies uji and uj2 satisfying uji + 0)2 = cop, up being the frequency of 
the pump beam, the present method can be used to generate squeezing at frequencies 
where no media for degenerate downconversion are available 

Phase-space approach has been also used to analyze CV teleportation as a 
conditional generalized double homodyning on a TWB. Also in this case the use of 
Wigner functions represents a powerful tool to evaluate the degrading effects of finite 
amount of entanglement, losses along the transmission channel, and nonunit quantum 
efficiency at sender location. A bound on the value of quantum efficiency needed to 
demonstrate quantum teleportation has been derived. 
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Figure 1. Schematic diagram of conditional homodyne for remote squeezed state 
preparation. A TWB is produced by spontaneous downconversion in a NOPA, 
and then homodyne detection is performed on one of the two modes, say mode 1. 
Mode 2 is squeezed if the homodyne quantum efficiency is larger than 50%. 




Figure 2. Schematic diagram of teleportation scheme. In the sender area (S), one 
part of a TWB is mixed with a given reference state a (the state to be teleported) 
in a balanced beam splitter, and two orthogonal quadratures are measured on 
the outgoing beams by means of two homodyne detectors with local oscillators 
phase-shifted by it/2. After the measurement, in the receiver area (R), the other 
part of the TWB is displaced by an amount —a = —x — iy that depends on the 
outcome of the measurements itself. The overall state, averaged over the possible 
outcomes, is the teleported state. 



